Nonlocal Effects on the Magnetic Penetration Depth in d-wave Superconductors 
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We show that, under certain conditions, the low temperature behavior ol the magnetic penetration 
depth A(T) ol a pure d-wave superconductor is determined by nonlocal electrodynamics and, contrary 
to the general beliel, the deviation AA(T) = A(T) — A(0) is proportional to T 2 and not T. We 
predict that the AA(T) oc T 2 dependence, due to nonlocality, should be observable experimentally 
in nominally clean high-T c superconductors below a crossover temperature T* — (£ /A ) A ~ IK. 
Possible complications due to impurities, surface quality and crystal axes orientation are discussed. 
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There is a significant amount of experimental evidence 
that the pairing state in the cuprate high temperature 
superconductors (HTSC) is unconventional, most proba- 
bly of d x 2_ y 2 symmetry jfl-pl. One generic feature of a 
layered HTSC with any unconventional order parameter 
(OP) compatible with the underlying crystal symmetry 
is that the OP exhibits line nodes (point nodes in 2D) on 
the Fermi surface (FS) and, therefore, gapless quasipar- 
ticle excitations in the corresponding energy spectrum. 
These low-lying excitations dominate the low tempera- 
ture thermodynamics and transport properties of these 
materials, and it is expected that the temperature de- 
pendence of the different thermodynamic quantities and 
transport coefficients will follow a power law rather than 
the conventional exponential behavior 8. Direct exper- 
imental evidence for the existence of zeros of the gap 
function on the FS in HTSC has been found by angle 
resolved photo emission spectroscopy (ARPES) |g,|). 

In particular, the low temperature behavior of the 
Meissner penetration depth X(T) is frequently regarded 
as an important probe of the morphology of the magni- 
tude of the anisotropic OP, A (p), in the cuprates. In con- 
ventional s-wave superconductors, the deviation AA(T) 
of A(T) from its zero temperature value A(0) exhibits 
activated behavior, i.e., AA(T) oc exp(— A/T) (through- 
out this paper we use units in which kg = fi = 1), re- 
flecting the existence of the isotropic BCS energy gap 
A at the FS. In contrast, in a pure d-wave supercon- 
ductor, or any other unconventional superconductor with 
nodes in the gap, the London (local) penetration depth 
varies linearly with the temperature, i.e., AA(T) oc T. 
Recently, by employing different high precision measure- 
ment techniques, such a linear T-dependence of the in- 
plane AA a (,(T) penetration depth (here the subscript 
refers to the axes along which the screening currents flow) 
has been observed experimentally in the Meissner state of 
several HTSC systems, such as: high quality single crys- 
tals of YBa 2 Cu 3 7 -5 (YBCO) @-| and Bi 2 Sr 2 CaCu 2 8 
(BSCCO) |lO|~|T^], magnetically aligned powders of crys- 
talline HgBa 2 Ca 2 Cu 3 8+5 JL3| and high quality YBCO 
thin films p4-TTq]. However, below a certain sample de- 



pendent temperature T* the linear T-dependence of the 
penetration depth in HTSC crosses over into a higher 
power law, most probably T 2 . In the d-wave scenario 
of HTSC, the origin of the A(T) oc T 2 dependence has 
been explained by the presence of non-magnetic impuri- 
ties which scatter in the unitary limit fll7| ] . In this strong 
scattering limit a small amount of impurities can induce 
a finite residual density of states at the Fermi level which 
is sufficient to change the temperature dependence of the 
penetration depth from T to T 2 without lowering signifi- 
cantly the transition temperature. A direct experimental 
confirmation of such a crossover between pure and impu- 
rity dominated regimes was reported by Bonn et al. [ fl8[ . 
The purpose of this Letter is to show that at very low 
temperatures nonlocality may play an important role in 
the electromagnetic response of a d-wave superconduc- 
tor (or any other unconventional superconductor with 
nodes in the gap), leading to a AA(T) oc T 2 depen- 
dence even in the clean limit. Thus, besides impuri- 
ties, nonlocality represents a second mechanism which 
leads to a T 2 dependence of the penetration depth suf- 
ficiently close to T = OK. To the best of our knowl- 
edge, all theoretical calculations and interpretations of 
the experimental measurements of the penetration depth 
in HTSC performed so far assume the validity of lo- 
cal electrodynamics (London limit) p9| . At first sight 
this is reasonable, since the zero temperature London 
penetration depth A = ^Jmc 2 /Aime 2 is much larger 
than the corresponding coherence length £ in these ma- 
terials. (In contrast to the penetration depth, which 
can be measured more or less directly, the coherence 
length £ cannot be determined experimentally and, in 
fact, is estimated in terms of the maximum value of the 
anisotropic gap function A = max{A(p)} by using the 
usual BCS expression £ = vp/ivA .) However, in the 
case of a clean, anisotropic superconductor it is more ap- 
propriate to introduce an anisotropic coherence length 
£ (p) = vp/ \A(p)\. If the anisotropic OP has nodes 
on the FS, it is clear that sufficiently close to the nodes 
£ (p) /A = (£ /A ) A / I A (p)| > 1 holds and, therefore, 
the contribution of these regions of the FS to the pene- 



tration depth X(T) must be determined by using nonlo- 
cal electrodynamics. The large value of A /£ guarantees 
that the applicability of local electrodynamics is violated 
only on a very small fraction, of order a = £ /A , of 
the FS. Since the whole FS contributes to the zero tem- 
perature penetration depth A(0), one expects no signifi- 
cant corrections to this quantity due to nonlocal effects. 
On the other hand, the low temperature dependence of 
A(T) must be dominated by nonlocal effects because this 
dependence is determined by a small region of the FS 
which is concentrated around the nodes of the OP. The 
crossover temperature below which nonlocal effects are 
important is given by T* — a A . Indeed, since the 
range of p values corresponding to the thermally excited 
quasiparticles at a given temperature T is determined 
by the condition |A(p)| < T, for T < T* one obtains 
f (p) /A = a A / |A (p)\ > T*/T > 1. For T » T* the 
local limit is applicable. As a typical example consider 
a YBCO single crystal with A sw 250K, £ sw 14A and 
A w 1400A; this yields a « I0~ 2 and T* w 2.5K. 

To demonstrate the effect of nonlocal electrodynamics 
on X(T) let us consider the case when a weak, uniform 
and static magnetic field H = V x A is applied along 
the c-axis of a semi-infinite HTSC with a plane surface 
which is perpendicular to the 6-axis. For this particular 
geometry, both the vector potential A and the screening 
supercurrent density j are oriented parallel to the a-axis, 
while the direction of penetration is along the 6-axis. We 
model the HTSC as a quasi two-dimensional rf-wave su- 
perconductor in which the motion of the electrons is con- 
fined for our purposes to the Cu-0 planes. In principle, to 
calculate A(T), one must first solve self-consistently the 
relevant Maxwell equation V x V x A = (Att/c) j together 
with the equation which relates j to A, subject to some 
properly chosen boundary conditions. In a weak mag- 
netic field (Meissner state) linear response theory yields 
for our geometry j(y) — — Jdy' K (y, y') A (y'), where y 
is the coordinate along the 6-axis (y = gives the position 
of the boundary) and the nonlocal electromagnetic re- 
sponse kernel K (y, y') must be calculated by using some 
microscopic theory. Once H(y) is determined, the pen- 
etration depth can be calculated according to the stan- 
dard definition (valid for a semi-infinite superconductor 
with plane boundary) A = iJ(0) _1 L dyH(y). Further- 
more, we assume that the boundary reflects the electrons 
either specularly or diffusively. In both these limiting 
cases X(T) can be expressed in terms of the Fourier trans- 
form of the bulk response kernel K(q; T). For a specular 
boundary one has [p0[ 



^specv^ ) 
Xo 



dq 



q 2 + K(q;T) 



while for a diffuse boundary [EG 
Adiff(T) 



A„ 



dq\n l + K(q;T)/(t 
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where the dimensionless quantities q and K are given by 
q = qX and K = [AirX^/c) K, respectively. 

For a weak-coupling, anisotropic superconductor the 
nonlocal bulk response kernel is similar to the corre- 
sponding expression for a conventional s-wave supercon- 
ductor |2 ill and can be written as 



K(q-T) = 2ttT Y, 




(3) 



where Lu n are fcrmionic Matsubara frequencies, A p = 
A (p), pit is the projection of p on the boundary, (. . .) 
means averaging over the circular 2D Fermi surface, and 
a = (qvp/2) qp. Here q is a unit vector perpendicular 
to the boundary and it gives the direction in which the 
penetration of the magnetic field takes place. Note that 
in a different geometry where the boundary is parallel to 
the a-b plane (H parallel to the boundary) , the direction 
of penetration q would be along the c-axis, i.e., perpen- 
dicular to p, yielding a = 0. Thus, we may conclude that 
the effect of nonlocal electrodynamics on X a b(T) is rele- 
vant only when H is parallel to the c-axis. Furthermore, 
at sufficiently low temperatures, the OP in Eq. (0) can 
be approximated with its limiting expression close to the 
nodes, i.e., A p = A <I> (p) w A o 4>'(0)<^, where <p is the 
angular deviation of p from the given node direction in 
the basal plane. In the case of a model d-wave OP with 
$ (p) =v\~v\ onc has $'(0) = 2. 

Let us calculate first the nonlocal correction to the zero 
temperature penetration depth A(0). For T = the fre- 
quency sum in (0) goes into an integral which can be 
evaluated exactly with the result 



K(q;0) 



2P|| 



sinh 1 (a/Ap) 



(a/A p ) yjl + (a/Apf 



(4) 



The average over the FS in (Q) can be evaluated ana- 
lytically for both London (local) and Pippard (extreme 
nonlocal) limits. In general, for arbitrary q values, 
K (q; 0) must be calculated numerically. In the Lon- 
don limit, when a q = q£ -C 1, one obtains K (q;0) = 
1 — (7r 2 v / 2/I6) ct q, while in the Pippard limit, when 
a q S> 1, one has K (q; 0) = (2/3) In (a q) / (a q) . Note 
that in both limiting cases the response kernel for a d- 
wave superconductor decreases with q more rapidly than 
for a conventional s-wave superconductor p0[ . Now the 
correction to the zero temperature penetration depth 
due to nonlocality can be obtained from Eqs. (|l[0) for 
both specular and diffuse boundaries. The results are: 



AA S p CC (0)/A o — A spoc (0)/A o 



Try/2ao/16, and 



AA d iff(0)/AA spcc (0) = In (o^ 2 ) /2 « 4.6. Thus for both 
type of boundaries, due to the very small value of a , the 
nonlocal correction to A(0) is less than 1% and therefore 
it can be obviously neglected, especially because this cor- 
rection is situated within the experimental errors of the 
most accurate measurements of the absolute value of the 
penetration depth. 

We turn now to calculate AA(T). At low temperatures, 
SK (q; T) = K (q; T) — K (q; 0) represents a small cor- 
rection to the zero temperature response kernel K (q; 0). 
Therefore, by using our previous result A(0) ~ A , from 
Eqs. (|Bj|) one obtains 



and 



AA spcc (T) 
A 



AA diff (r) 

A 



2 - ry**®? 

Wo {q 2 + lf 



1 f°° r ,-SK(q;T) 

dq — - 2 



q 2 + l 



(5) 



(6) 



Furthermore, a convenient expression for SK can be 
obtained by evaluating the Matsubara sum in Eq. (||) by 
means of complex contour integration 



SK(q;T) 



f(w) duj 



2pf, Re 



A 2 (A 2 
p \ v 



(7) 



Note that in the a — > limit one recovers the familiar 
local limit expression for SK |17|| . Due to the presence 
of the Fermi function f(u>) in (CT the main contribution 
to the frequency integral comes from the interval u> ;$ T. 
Therefore, in the average over the FS the relevant re- 
gions are determined by |A p | <uj^T and are obviously 
located around the nodes of the OP. By using the expres- 
sion for the OP close to a node one arrives, after some 
straightforward algebra, at the following result 



SK(q;T) = SK(0;T)F 



(8) 



where t = T/T*, SK(0;T) = -21n2T/A is the well 
known local limit expression of SK for a d-wave super- 
conductor yj], the expression 



F(z) = 1- 



ZLVf , 
1 /• — 2 

^J 



dxf{x)Jl-±(Z) (9) 



is a universal function, and f(x) = (e x + 1) .It is 
remarkable that, within the above mentioned approxi- 
mations, the kernel SK depends only on the ratio q/t 
and not separately on its two arguments. In order to 



make further analytical progress it is desirable to ap- 
proximate F(z) by a simpler expression. A careful quan- 
titative analysis of Eq. (||) motivates the following rea- 
sonable approximation: F(z) « 1 — c\z, for z < 2, and 
F(z) w c /z 2 , for z > 2, where c = 6C(3)/7r 2 ln2 « 1.05, 
and ci = (1 — c /4) /2 w 0.37. Note that ci is somewhat 
smaller than the absolute value of the slope of F(z) at 
the origin, i.e., |F'(0)| = tt 2 \/2/32 ln2 « 0.63. The tem- 
perature dependence of the penetration depth can now 
be calculated by inserting (||) in Eqs. (|^-||). 

For t> 1 (i.e., T ~^> T*) and for a specularly reflect- 
ing boundary one obtains AA S p ec (? 1 )/A = In 2 (T/A ) — 
(-71-^/2/16) a + 0(l/t). Here, the leading term is the 
well known linear in T local expression for AA(T) for a 
d-wave superconductor, i.e., AAl(T) = In 2 (T/A ) A . 
The second, small negative constant term in the ex- 
pression of AA spoc (T) is due to nonlocality and shows 
clearly that the linear T-dependence cannot extend all 
the way down to T — OK; it must cross over to a 
higher power law at some T ~ T* . In the case of a 
diffuse boundary one obtains a similar result, namely 
AA diff (T) = AA L (T) - (?r>/2/16) a X kit + 0(l/t). 
Note that the magnitude of the nonlocal correction to 
the local penetration depth is larger than in the case of 
the specular boundary by a factor of Int. 

In the opposite limit t <C 1 (i.e., T <^T*) one obtains 
for a specular boundary AA spoc (T) = (3 A\ L (T)T/T* ex 
T 2 , where /? = 8 (1 - c x + c /4) /tt « 2.2. Thus, due to 
nonlocal electrodynamics, for T <^i T* the temperature 
dependence of a pwre d-wave superconductor is propor- 
tional to T 2 and not T, regardless of how small is a = 
^ /A . This conclusion is one of the main results of the 
present paper. It should be noted that the above value 
for the coefficient /3 is just an approximation; a more ac- 
curate value of (3 can be obtained by approximating F(z) 
by a polynomial of degree N > 1 for z < z , and by its 
large z asymptotic form for z > z Q , where z is a con- 
veniently chosen value. By reevaluating the integrals in 
AA(T) one obtains again the leading term proportional to 
T 2 but with a slightly different numerical value for j3. A 
similar calculation in the case of a diffuse boundary yields 
AA diff (T) w AA spoc (T)/2. Note that in both limiting 
cases the inequalities AAdig(r) < AA spcc (T) < AAl(T) 
hold and imply that the deviation from the correspond- 
ing local result is larger for a diffuse boundary than for 
a specular one (see also Fig. |l|). 

For arbitrary temperatures AA(T) must be calculated 
numerically by employing the exact expression (0) for 
the function F(z). In Fig. |a the ratio AA(T)/AA L (T) 
is plotted, for both specular and diffuse boundaries, as 
a function of the reduced temperature t. The deviation 
from the standard result obtained in the local limit is ev- 
ident. The clear linear dependence in the vicinity of the 
origin indicates a quadratic T-dependence of AA(T). For 
t 3> 1, AA(T) approaches asymptotically its local limit 
(minus a small constant correction of order a )- 
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FIG. 1. Plot of AA(T) [in units (a) AX L (T), and (b) a \ , 
respectively] vs. t = T/T* for both specular (solid line) and 
diffuse (long-dashed line) boundary. For comparison, the lo- 
cal limit result is also shown (dashed line). 



high-T c superconductor. At temperatures T <C T* ~ IK 
the penetration depth X(T) has a quadratic temperature 
dependence, while above the crossover temperature T* , 
but still well below T c , A(T) has the well known linear T- 
dependence. Thus, nonlocality represents a second pos- 
sible mechanism, beside strongly scattering impurities, 
which may account for the experimentally observed de- 
viation from the linear T-dependence of the penetration 
depth at the lowest measured temperatures in nominally 
clean HTSC. A simple experiment to probe the viability 
of this mechanism has been proposed. 
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helpful discussions. One of us (A.J.L.) thanks G. E. 
Volovik for a conversation in which the possibility of this 
effect first emerged. This work was supported by the Na- 
tional Science Foundation (DMR 91-20000) through the 
Science and Technology Center for Superconductivity. 



Note that the deviation of AA(T) from the correspond- 
ing local expression is much more pronounced for a diffuse 
boundary then for a specular one. The same AA(T), this 
time in units of a \ , is shown as a function of t = T/T* 
in Fig. |ljb. The deviation from linearity becomes visible 
around t = 1 (t = 2) for the specular (diffuse) boundary. 
For the numerical example considered above for a clean 
YBCO single crystal (a X — £ « 14A) one finds that 
the deviation from AA(T) oc T takes place somewhere 
between 2 to 5 K, depending on the surface quality of 
the crystal. Such a crossover is seen experimentally in 
nominally clean YBCO crystals (22|; in the existing liter- 
ature it has been attributed to impurities which scatter 
in the unitary limit |17|. 

In principle, there is a simple experimental test to de- 
termine whether this crossover in AA(T) is due to non- 
local electrodynamics or to impurities. The main idea is 
to estimate experimentally the crossover temperature in 
AAab(T') for the same nominally clean HTSC for two dif- 
ferent magnetic field orientations: (i) H parallel to the 
c-axis, and (ii) H parallel to the a-b plane. As we have 
already mentioned, nonlocality is expected to be relevant 
only when the applied magnetic field is oriented parallel 
to the c-axis (so that the penetration direction lies in 
the a-b plane), while the effect of impurities should not 
depend on the orientation of the field. Thus, if T* is 
noticeably smaller in case (ii) than in case (i) one may 
conclude that the observed effect is mainly due to non- 
local electrodynamics and not to impurities. Otherwise, 
the conclusion is that nonlocal effects are in fact com- 
pletely masked by impurities. 

In conclusion, we have shown that nonlocal electro- 
dynamics dominate the low temperature behavior of the 
in-plane magnetic penetration depth of a clean d-wave 



[1] D. J. Scalapino, Phys. Rep. 250, 329 (1995). 
[2] D. J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 
[3] J. F. Annett, N. Goldenfeld, and A. J. Leggett, in 
Physical Properties of High Temperature Superconduc- 
tors, edited by D. M. Ginsberg (World Scientific, Sin- 
gapore, 1996), Vol. 5, Chap. 6, pp. 375-461. 
[4] M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
[5] H. Ding et al, Phys. Rev. B 50, 1333 (1994). 
[6] H. Ding et al, Phys. Rev. Lett. 74, 2784 (1995). 
[7] W. N. Hardy et al, Phys. Rev. Lett. 70, 3999 (1993). 
[8] K. Zhang et al, Phys. Rev. Lett. 73, 2484 (1994). 
[9] J. Mao et al, Phys. Rev. B 51, 3316 (1995). 
[10] T. Jacobs et al, Phys. Rev. Lett. 75, 4516 (1995). 
[11] S. F. Lee et al., Phys. Rev. Lett. 77, 735 (1996). 
[12] O. Waldmann et al, Phys. Rev. B 53, 11825 (1996). 
[13] C. Panagopoulos et al, Phys. Rev. B 53, R2999 (1996). 
[14] O. M. Froehlich et al, Phys. Rev. B 50, 13894 (1994). 
[15] O. M. Froehlich et al, Euro Phys. Lett. 36, 467 (1996). 
[16] A. Fuchs, W. Prusseit, P. Berberich, and H. Kinder, 

Phys. Rev. B 53, 14745 (1996). 
[17] J. P. Hirschfeld and N. Goldenfeld, Phys. Rev. B 48, 4219 

(1993). 
[18] D. A. Bonn et al, Phys. Rev. B 50, 4051 (1994). 
[19] Nonlocal electrodynamics has been applied recently to 
study the optical response of ultra-clean d-wave HTSC 
by J. A. Sauls (unpublished) and D. Xu, Ph.D. thesis, 
Northwestern University, 1995, (unpublished). 
[20] M. Tinkham, Introduction to Superconductivity (McGraw 

Hill, New York, 1974). 
[21] A. A. Abrikosov, L. P. Gor'kov, and I. E. Dzyaloshin- 
ski, Methods of quantum field theory in statistical physics 
(Prentice-Hall, Englewood Cliffs, N.J., 1963). 
[22] D. A. Bonn and W. N. Hardy, in Physical Properties 
of High Temperature Superconductors, edited by D. M. 
Ginsberg (World Scientific, Singapore, 1996), Vol. 5, 
Chap. 2, pp. 7-98. 



